Establishment and Characterization of Immortalized Human Dermal Papilla Cells Expressing Human Papillomavirus 16 E6/E7

Primary human dermal papilla cells (HDPCs) are often preferred in studies on hair growth and regeneration. However, primary HDPCs are limited by their reduced proliferative capacity, decreased hair induction potential, and extended doubling times at higher passages. To overcome these limitations, pTARGET vectors containing human papillomavirus16 (HPV16) E6/E7 oncogenes were transfected into HDPCs and selected using G-148 to generate immortalized cells here. HPV16 E6/E7 oncogenes were efficiently transfected into primary HDPCs. Immortalized HDPC showed higher proliferative activity than primary HDPC, confirming an increased proliferation rate. Expression of p53 and pRb proteins was downregulated by E6 and E7, respectively. E6/E7 expressing HDPC cells revealed that cyclin-dependent kinase (CDK) inhibitor p21 expression was decreased, while cell cycle-related genes and proteins (CDK2 and cyclin E) and E2F family genes were upregulated. Immortalized HDPCs maintained their responsiveness to Wnt/β-catenin pathway and hair follicle formation capability, as indicated by their aggregative properties and stemness. E6/E7 immortalized HDPCs may facilitate in vitro hair growth and regeneration studies.


Introduction
Hair dermal papilla cells (HDPCs) derived from the mesenchyme reside at the base of hair follicles.HPDCs play a crucial role in hair growth and development through interactions with the epithelial cells in hair follicles [1,2].The functions of HDPCs should be preserved to support continuous hair growth and development.However, cultured primary HDPCs lose hair inductive and proliferative capacity at higher passages [3].Thus, establishing immortalized HDPCs that retain the properties of primary HDPCs is necessary for hair growth and regenerative research.
Transfection of cell lines with HPV16 genes, as in HPV infection of normal epithelial cells or transduction with lentiviral vectors harboring HPV E6/E7 genes, is an established method for cell line immortalization [4].HPV16 E6/E7 oncogenes have been used to immortalize various cell lines, including human ameloblastoma, human oral epithelial, and bovine mammary epithelial cells [5][6][7].HPV is a small, non-enveloped, double-stranded DNA virus that encodes early (E1, E2, E4, E5, E6, and E7) and late (L1 and L2) proteins [8].Of these proteins, E6 and E7 oncoproteins from high-risk HPV promote cellular proliferation by deactivating p53 and retinoblastoma (Rb), respectively [9,10].E6 plays a central role in cell line immortalization by forming a ternary complex with the ubiquitin ligase E6 associated protein (E6AP) and p53 to trigger ubiquitination and subsequent degradation of p53 [11].E7 is primarily located in the nucleus, where it binds and phosphorylates pRb, resulting in the release of E2F [12].Immortalized HDPCs can be established by inactivating these major tumor suppressor proteins, p53 and pRb, through E6/E7 transfection [13].
Here, we transfected HPV16 E6/E7 oncogenes into early-passage primary HDPCs to generate immortalized HDPCs.We also performed an in vitro characterization of the immortalized HDPCs in terms of their proliferative capacity, Wnt/-catenin signaling pathway activity, and hair follicle formation capability.

Cell Culture
Primary human follicular dermal papilla cells (HDPCs) were obtained from PromoCell GmbH (Germany).The cells were cultured in hair follicle dermal papilla cell growth medium supplemented with 40 μl/ml of fetal calf Primary human dermal papilla cells (HDPCs) are often preferred in studies on hair growth and regeneration.However, primary HDPCs are limited by their reduced proliferative capacity, decreased hair induction potential, and extended doubling times at higher passages.To overcome these limitations, pTARGET vectors containing human papillomavirus16 (HPV16) E6/E7 oncogenes were transfected into HDPCs and selected using G-148 to generate immortalized cells here.HPV16 E6/E7 oncogenes were efficiently transfected into primary HDPCs.Immortalized HDPC showed higher proliferative activity than primary HDPC, confirming an increased proliferation rate.Expression of p53 and pRb proteins was downregulated by E6 and E7, respectively.E6/E7 expressing HDPC cells revealed that cyclin-dependent kinase (CDK) inhibitor p21 expression was decreased, while cell cycle-related genes and proteins (CDK2 and cyclin E) and E2F family genes were upregulated.Immortalized HDPCs maintained their responsiveness to Wnt/β-catenin pathway and hair follicle formation capability, as indicated by their aggregative properties and stemness.E6/E7 immortalized HDPCs may facilitate in vitro hair growth and regeneration studies.

RGB Profiling of Immunofluorescence Images
Images were prepared to include a legend with a scaling bar to represent actual cell sizes.The scale was adjusted using the ImageJ software.Briefly, a representative cell was selected and a line was drawn at the location that passed through the nucleus.The fluorescence intensity along the drawn line was measured using the red-greenblue (RGB) profiler of ImageJ software.This generates a single plot that shows the profile of an RGB image in three different colors (red, green, and blue).

Evaluation of 3D Spheroid Construction Ability
To assess whether HDPCs induce the 3D spheroid, primary and immortalized HDPCs (1.5 × 10 5 cells/ml) were seeded in a 96-well round-bottom cell floater plate.The cells were incubated for 48 h to generate spheroids.Images of the 3D spheroids were obtained using a phase-contrast microscope.The average value of each spheroid was estimated by measuring the diameter five times using the ImageJ software.

Statistical Analysis
Data were presented as the mean standard deviation (SD) of three independent experiments.GraphPad Prism 9 software (GraphPad Software Inc., USA) was used for the statistical analysis.Unpaired t-tests were used.Statistical significance was defined as follows: *p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001.

Detection of HPV16 E6/E7 in HDPCs
The primary HDPC cells were stably transfected with either E6/E7 from the integrated HPV16 genome or a mock vector, and the effect of E6/E7 on cell cycle progression was investigated.RT-PCR was used to quantify the expression of E6 and E7 in the respective transfectants (Fig. 1A and 1B).HPV16 E6 and E7 mRNA levels were significantly higher in immortalized HDPCs compared to primary HDPCs.Most E6 transcripts were full-length E6, and the spliced form E6*I was detected at low levels.

The Effect of E6/E7 on Proliferation and Morphology in HDPCs
The proliferation rates of primary and immortalized HDPCs were determined using the MTS assay.The rate was measured every 24 h for 72 h.The proliferation rates of primary and immortalized HDPCs were 149.7% and 255.3%, respectively.Hence, HPV16 E6/E7 immortalized HDPCs proliferated at a 1.71-fold higher rate than primary HDPCs under the same conditions (Fig. 2A).Both primary and immortalized HDPCs showed uniform cell sizes and did not show any differences from each other (Fig. 2B).These results indicate that E6/E7 oncoproteins increase the proliferation rate without affecting cell characteristics.

The Effect of E6 on Expression of p53 and Cell Cycle Related Factors
First, p53 levels were quantified in primary and immortalized HDPCs.Gene expression levels of TP53 were not altered (Fig. 3A), whereas expression of p53 protein was downregulated (Fig. 3B).Downregulation of p53 protein expression in immortalized HDPCs, especially in the nucleus, was also confirmed using immunofluorescence staining (Fig. 3C).We then assessed the effect of E6/E7 on the cyclin-dependent kinase (CDK) inhibitor p21.p21 protein expression was downregulated in immortalized HDPCs (Fig. 3D).CDK2 and CCNE1 mRNA levels were upregulated in immortalized HDPC (Fig. 3E).Protein expression levels of CDK2 and cyclin E were also increased in E6/E7 immortalized HDPCs (Fig. 3A and 3B).These results indicate that p53 is degraded through the proteasome via E6, leading to the decrease in p21 expression, and increase in CDK and cyclin E.

The Effect of E7 on Expression Levels of pRb and Transcription Factor E2F
Next, the direct effects of E7 on pRb and E2F expression in HDPCs were evaluated.Immunofluorescence staining revealed that pRb expression was downregulated in immortalized HDPCs, especially in the nucleus (Fig. 4A).The reduction in pRb expression and increase in p-pRb levels in immortalized HDPCs were also confirmed by immunoblotting assays (Fig. 4B and 4C).In addition, mRNA levels of E2F factors, including E2F1, E2F2, and E2F3, were increased in immortalized HDPCs (Fig. 4D).These results indicate that E7 hyperphosphorylates pRb, leading to the release of E2F which promotes transcription of genes involved in G1-S transition.

Responsiveness to Wnt/β-Catenin Signaling and Hair Follicle Induction in E6/E7 Immortalized HDPCs
Primary and immortalized HDPCs were treated with DHT (1 μM) to assess their Wnt/-catenin pathway responsiveness.DHT treatment for 2 h inhibited the nuclear translocation of β-catenin in both HDPC types   (Fig. 5A).By RGB line profiling, it was also confirmed that β-catenin intensity in the nucleus was low (Fig. 5B).These findings indicate that immortalized HDPCs retained their responsiveness to the Wnt/-catenin signaling pathway.
Since the induction of hair follicles is related to the aggregation of HDPCs, we also investigated their 3D spheroid formation capability of the transfected cells.Both primary and immortalized HDPCs formed round and compact 3D spheroids with a size of approximately 420 μm (Fig. 5C and 5D).To confirm whether E6/E7 affects the stemness of HDPCs, epithelial-mesenchymal transition (EMT) and self-renewal markers were quantified.EMT markers were well expressed even in the immortalized HDPCs (Fig. 5E).Immortalized HDPCs showed a slight decrease in SOX2 expression, whereas there were no noticeable differences in NANOG expression (Fig. 5F).These results indicated that both the aggregation ability and stem cell-like characteristics were preserved in immortalized HDPCs.
HPV E6*I binds p53 to induce proteasomal degradation of p53, which may result in cell immortalization [18].p53 regulates cellular responses including DNA repair, cell cycle arrest, apoptosis induction, and metabolism [19][20][21][22].TP53 expression did not differ between primary and immortalized HDPCs (Fig. 3A).However, p53 protein levels were lower in immortalized HDPCs compared to primary HDPCs (Fig. 3B-3C).These results indicate that p53 underwent increased proteasomal degradation through binding to the E6 protein.The CDK inhibitor p21 is the transcriptional target of p53.When activated, p53 binds the p21 promoter to increase transcription [23,24].As expected, p21 protein expression was decreased due to the downregulated p53 (Fig. 3D).p21 regulates cell cycle progression by inhibiting CDK proteins, which form complexes with cyclin to control various cell cycle stages [25].The G1-S phase transition is associated with CDK2, which is activated by cyclin E1 and E2 [26,27].CDK2 and CCNE1 mRNA levels and CDK2 and cyclin E protein levels were significantly upregulated in immortalized HDPCs (Fig. 3E-3G).These upregulations promote G1-S transition and induce cell proliferation.
E7 in high-risk HPVs has a higher affinity for pRb than E7 in low-risk HPVs, leading to a reduction in pRb level [28].E7 interacts with retinoblastoma (Rb) family proteins, such as pRb, p107, and p130 [29].pRb expression was decreased in immortalized HDPCs (Fig. 4A-4B).However, p-pRb expression was increased compared to primary HDPCs (Fig. 4C), which was due to the increase in CDK and cyclin E as shown in Fig. 3E-3F.Upon formation in G1, cyclin E-CDK2 complexes translocate to the nucleus and phosphorylate pRb, effectively inactivating it [30].Free E2F regulates the expression of genes, including cyclin, CDK, and E2F themselves, which regulate the G1-S transition [26,31].Only E2F1-3 among the E2F family bind pRb, and its overexpression efficiently induces G1-S transition [28,32].We measured E2F1, E2F2, and E2F3 mRNA levels in the two types of HDPCs and found significant upregulation in immortalized HDPCs (Fig. 4D).This upregulation was expected to increase the expression of genes critical for G1-S transition.
As a potent androgen, DHT plays a key role in regulating the hair cycle.Hair cycling relies on crosstalk between the androgen and Wnt/β-catenin pathways [33].Wnt/β-catenin pathway is active in the dermal papilla of growing hair follicles during the anagen phase, and maintaining the activity of this signaling pathway is essential for preserving cultured HDPCs in a state that promote hair growth [34].β-catenin activity was inhibited by DHT treatment (Fig. 5A and 5B), confirming that the response to Wnt/β-catenin signaling pathway sustained functional even in immortalized HDPCs.
HDPCs are multipotent stem cells and their aggregative properties and stemness are associated with hair follicle formation [35,36].3D spheroids efficiently induce HDPC that can form hair follicles in the skin [37].Previous studies demonstrated that 3D spheroids of HDPCs form de novo hair follicles on mice skin [38].Models with cocultures incorporating both dermal papilla cells and keratinocytes have also been developed, providing evidence that their co-cultivation effectively preserves the in vivo properties of dermal papilla cells [39,40].The human folliculoid microsphere, consisting of HDPC and outer root sheath keratinocytes, has also been reported.This model demonstrated close interaction between the two types of cells and expressed versican, which is essential for hair follicle formation [41].In addition, maintenance of the stem cell-like characteristics is crucial to maintain HDPC functionality.EMT and self-renewal factors are markers of stemness in HDPCs [35].Transcription factors such as Twist, Snail, and Slug are found in cells undergoing EMT and play a significant role in preserving stem cell function [42,43].Self-renewal, a key characteristic of stem cells, is facilitated by SOX2 and NANOG [44].SOX2 is a well-known HDPCs marker [45].Thus, the aggregative behavior, EMT markers (TWIST, SNAIL, SLUG, and ZEB1), and self-renewal markers (SOX2 and NANOG) were evaluated in both cell types.E6/E7 transfection did not affect the aggregation or stemness of HDPCs (Fig. 5C-5F).This suggests that immortalized HDPCs retained their ability to form hair follicles while maintaining their characteristic features.

Conclusion
Here, we established and characterized immortalized HPDC cells by transfecting them with the HPV 16 E6/E7 oncogenes.Upon transfection, proliferation rate of cells was promoted through regulating of p53 and pRb activity.The cells also exhibited preserved responsiveness to the Wnt/β-catenin pathway, and retained their ability to form hair follicles, as demonstrated by their aggregating properties and stem cell characteristics.Immortalized HDPCs may be useful for in vitro research on mechanisms of hair growth and regeneration.

Fig. 3 .
Fig. 3. Expression levels of p53 and cell cycle related factors in HDPCs.(A, B) The TP53 gene and p53 protein expression level was determined.The relative band intensities were evaluated by ImageJ software.Data are presented as the mean SD (n = 3).***p < 0.001 by t-test.(C) p53 was labeled with Cy3 (red) in cells, and nuclei were labeled with DAPI (blue).Scale bar = 125 μm.(D) The protein level of CDK inhibitor p21 was determined by immunoblotting assay.(E-G) The gene (CDK2 and CCNE1) and protein (CDK2 and cyclin E) expression levels were determined.The relative band intensities were evaluated by ImageJ software.Data are presented as the mean SD (n = 3).*p < 0.05, **p < 0.01 by t-test.

Fig. 4 .
Fig. 4. Expression levels of pRb protein and E2F genes in HDPCs.(A) pRb was labeled with FITC (green) in cells, and nuclei were labeled with DAPI (blue).Scale bar = 125 μm.(B, C) The pRb and p-pRb protein level was determined by immunoblotting assay.The relative band intensities were evaluated by ImageJ software.Data are presented as the mean SD. **p < 0.01 by t-test.(D) The gene expression levels of E2F1-3 were determined.The relative band intensities were evaluated by ImageJ software.Data are reported as the mean SD (n = 3).*p < 0.05, **p < 0.01, ***p < 0.001 by t-test.

Fig. 5 .
Fig. 5. Responses to Wnt/β-catenin signaling and hair follicle formation capability in HDPCs.(A) β-catenin was labeled with Cy3 (red) in cells, and nuclei were labeled with DAPI (blue).Scale bar = 125 μm.(B) Fluorescence intensity was profiled (white arrow sign on third line of (A)).(C, D) Representative image of 3D spheroid was confirmed (magnification at × 4).Scale bar = 500 μm.The size in diameter of spheroids were quantified by ImageJ software.(E, F) The gene expressions of EMT markers (TWIST, SNAIL, SLUG, and ZEB1) and self-renewal markers (SOX2 and NANOG) were estimated.The relative band intensities were evaluated by ImageJ software.Data are reported as the mean SD (n = 3).**p < 0.01 by t-test.